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Through exact numerical solutions we show Bose-Einstein condensation (BEC) for the one- 
dimensional (ID) bosons with repulsive short-range interactions at zero temperature by taking 
a particular large size limit. Following the Penrose-Onsager criterion of BEC, we define condensate 
fraction by the fraction of the largest eigenvalue of the one-particle reduced density matrix. We 
show the finite-size scaling behavior such that condensate fraction is given by a scaling function of 
one-variable: interaction parameter multiplied by a power of particle number. Condensate fraction 
is nonzero and constant for any large value of particle number or system size, if the interaction 
parameter is proportional to the negative power of particle number. Here the interaction parameter 
is defined by the coupling constant of the delta-function potentials devided by the density. With 
the scaling behavior we derive various themodynamic limits where condensate fraction is constant 
for any large system size; for instance, it is the case even in the system of a finite particle number. 
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The experimental realization of trapped atomic gases 
in one dimension has provided a new motivation for the 
study of strong correlations in fundamentalquantum me- 
chanical systems of interacting particles [l|-|3[. In one- 
dimensional (ID) systems quantum fluctuations play a 
key role and often give subtle and nontrivial effects. It is 
known that Bose-Einstein condenstation (BEC) occurs 
even for bosons with repulsive interactions due to the 
quantum statistical effect among identical particles [3]. 
In fact, the existence of BEC has been proven rigorously 
for interactiong bosons confined in dimensions greater 
than one [Hj. In ID case there is no BEC for bosons with 
repulsive interactions due to strong quantum fluctuations 
if we take the standard thermodynamic limit with fixed 
coupling constant On the other hand, if the coupling 
constant is very weak, we may expect that even the ID 
bosons with a large but finite number of particles undergo 
a quasi-condensation in which "a macroscopic number of 
particles occupy a single one-particle state" [3] . 

However, it has not been shown explicitly how such 
a quasi-condensation occurs in interacting bosons in one 
dimension. Furthermore, it is nontrivial to expect it for 
the ID Bose gas that is solvable by the Bethe ansatz. 
No pair of particles can have the same quasi-momentum 
in common for a Bethe-ansatz solution. Here we call 
the ID system of bosons interacting with repulsive delta- 
function potentials the ID Bose gas. For the impene- 
trable ID Bose gas where the coupling constant is taken 
to infinity, condensate fractions are analytically and nu- 
merically studied 0, while in the weak coupling case it 
is nontrivial to evalaute the fractions in the ID Bose gas. 

In the present Letter we show the finite-size scaling be- 
havior of condensate fraction in the ID Bose gas with re- 
pulsive interactions at zero temperature, and derive BEC 
by controlling the thermodynamic limit. We show that if 
the coupling constant decreases as a power of the system 
size, condensate fraction does not vanish and remains 



constant when we send the system size to infinity or to a 
very large value with fixed density. If condensate fraction 
is nonzero for a large number of particles, we call it BEC 
according to the Penrose-Onsager criterion. 

The scaling behavior of BEC in the ID Bose gas is 
fundamental when we specify the thermodynamic limit, 
where we send particle number N or system size L to in- 
finity or very large values. We define interaction param- 
eter 7 by 7 = c/n with coupling constant c in the delta- 
function potentials and density n = N/L. We show that 
if 7 is given by a negative power of iV, i.e. 7 = A/N^ , 
condensate fraction ng is nonzero and constant for any 
large value of L or N. We also show that exponent 77 
and amplitude A are independent of density n, and eval- 
uate them as functions of hq. Condensate fraction no is 
thus given by a scaling function of variable jN^. If the 
condensate fraction of a quantum state with large is 
nonzero in the ID Bose gas, we suggest that the classical 
mean-field approximation such as the Gross-Pitaevskii 
(CP) equation is valid for the state Furthermore, we 
show that the ID Bose gas of a finite particle number 
may have the same condensate fraction for any large L. 

Let us review the definition of BEC through the one- 
particle reduced density matrix for a quantum system 
[1, 0]. We assume that the number of particles is 
very large but finite. At zero temperature, the density 
matrix is given by p = |A)(A|, where |A) denotes the 
ground state of the quantum system. We define the one- 
particle reduced density matrix by the partial trace of the 
density matrix with respect to other degrees of freedom: 
pi = NtT23---NP- This matrix is positive definite and 
hence it is diagonalized as 



pi = Afo|«'o)(«'o| 



(1) 



Here we put eigenvalues Nj in descending order: A^o > 
A^i > A^2 > • • • > 0. The sum of all the eigenvalues is 
given by the number of particles: ^ Nj = N. Here we 
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recall tripi = N due to the normalization: tri23...7v/5 = 1- 
Let us denote by uq the ratio of the largest eigenvalue Nq 
to particle number N: 



no := Na/N. 



(2) 



The criterion of BEC due to Penrose and Onsager [9| is 
given as follows: If the largest eigenvalue Nq is of order 
N, i.e. the ratio no is nonzero and finite for large N, 
then we say that the system exhibits BEC, and we call 
no the condensate fraction. Here we also define fractions 
rij by 7ij = Nj/N for j = 1, 2, . . .. 

We now consider the Hamiltonian of the ID Bo se g as, 
which we call the Lieb-Liniger model (LL model) [10|: 



N ^2 

V — 



N 

j<k 



Xk)- 



(3) 



We assume the periodic boundary conditions of the sys- 
tem size L on the wavefunctions. We employ a system 
of units with 2m = h = 1, where m is the mass of the 
particle. We consider the repulsive interaction: c > 0, 
hereafter. 

In the thermodynamic limit the LL model is charac- 
terized by the single parameter 7 — c/n^ where n = N/L 
is the density of particle number N. We fix the particle- 
number density as n = 1 throughout the Letter, and 
change coupling constant c so that we have different val- 
ues of 7. 

In the LL model, the Bethe ansatz offers an exact 
eigenstate with an exact energy eigenvalue for a given 
set of quasi-momcnta fci, fc2, • ■ ■ , satisfying the Bethe 
ansatz equations (BAE) for j = 1, 2, . . . , A^: 



kjL 



27r/, 



N 

2 arctan 



ki — kf 



(4) 



Here Ij 's are integers for odd N and half-odd integers for 
even N. We call them the Bethe quantum numbers. The 
total momentum P and energy eigenvalue E are written 
in terms of the quasi-momenta as 



N 



2tt 



N 



E 



N 

i=i 



If we specify a set of Bethe quantum numbers Ii < ■ ■ ■ < 
In, BAE (jH) have a unique real solution fci < • ■ • < fcjv 
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In particular, the sequence of the Bethe quan- 
tum numbers of the ground state is given by Ij — 
-{N + l)/2 + j for integers j with I < j < N. The 
Bethe quantum numbers for low lying excitations are sys- 
tematically derived by putting holes or particles in the 
perfectly regular ground-state sequence. 

The matrix element of the one-particle reduced density 
matrix, pi(x,y) := (x\pi\y), for a quantum system is 



0.01 



100 



Iplh 0.999984 0.971538 0.693620 
2p2h 1.59454 x 10"^ 0.0280102 0.289056 
nsat 1.00000 0.999548 0.982676 

TABLE I: Fraction risat of the reduced density operator 
at the origin, pi (0,0), to the density n, evaluated by tak- 
ing the sum over a large number of eigenstates |^) with 
one particle and one hole (Iplh) or with two particles and 
two holes (2p2h) for iV = L = 50 (n = 1): riaat = 

Er'+Er')i(Mi^(o)iA)ivn. 



expressed as a correlation function in the ground state 
|A): 



(6) 



In the LL model we can numerically evaluate the corre- 
lation function by the form factor expansion. Inserting 
the complete system of eigenstates, |/i)(/i|, we have 



Piix^y) 



I(mI^(o)|a)P, (7) 



where denote the momentum eigenvalues of eigen- 
states 111). Each form factor in the sum ^ is expressed 
as a product of determinants by making use of the deter- 
minant formula for the norms of Bethe eigenstates 11211 
and that for the form factors of the field operator 13l4l5j : 



{fi\m\^) = (-1) 

=1 e=i 1 



JV(Af+l)/2+l 




k'. 



det U{k, k') 



^detG(fc) detG'(fc') 



(8) 



where 



the quasi-momcnta {fci,-- - ,fcAr} and 
i^Ar-i} Si'^^ eigenstates |A) and |/x), re- 
spectively. Here we have employed the abbreviated 
symbols kjj :~ kj — kf and k'^ ^ :~ k'j — kg. The matrix 
G{k) is the Gaudin matrix, whose (j, ^)th element 



(5) is G(fc),-f 



L + El=lKik,,m) 

j,£ = 1,2,- •• ,A^, where the kernel K{k) is defined by 
K{k) = 2c/(fc2 + c2). The matrix elements of the (A^-1) 
by (A^ — 1) matrix U {k, k') arc given by [l^ - [l5l | 



K{kjj) for 



ic) 



ic) 



{K{kj^,) - K{kN,l)) . (9) 



Numerically we calculate correlation function ([7]) by 
taking the sum over a large number of eigenstates with 
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FIG. 1: (Color online) Dependence of fractions Uj on coupling 
constant c. In the upper panel: condensate fraction no is 
plotted against coupling constant c for A*' = 4, 10, 20, 40, 
100, 200 and 400, from the top to the bottom, in red, green, 
blue, black, orange, purple and cyan lines, respectively. In the 
lower panel: condensate fraction tiq, fractions ni and 712 are 
shown against c from the top to the bottom in blue, red and 
green lines, respectively, for A'^ = 20. We recall n = N/ L = 1. 



one particle and one hole (Iplh) and those with two par- 
ticles and two holes (2p2h). In order to confirm the valid- 
ity of the restricted sum, we have estimated the ratio of 
the one-particle reduced density operator at the origin to 
density 71, pi(0, 0)/?!, through the form factor expansion 
([7|) for the excitations with Iplh or 2p2h. We denote it 
by risat- The estimates of risat are listed in Table 1. The 
graph of risat approaches 1 for small coupling constant c, 
while it is larger than 0.98 for any value of c in the case 
of TV = 50. 

For the LL model, the eigenfunctions of the one- 
particle reduced density matrix are given by plane waves 
for any nonzero and finite value of c. It is a consequence 
of the translational invariance of the Hamiltonian of the 
LL model. We thus have 



2N- 

.^_lcos[27rj(x-y)/L] . (10) 



The eigenvalues of the one-particle reduced density ma- 
trix, iVj, are expressed in terms of the form factor ex- 
pansion. We consider the sum over all the form factors 
between the ground state, |A), and such eigcnstates, |^), 
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FIG. 2: (Color online) Condensate fraction no as a function 
of l/N for c = 0.01. Here n = N/L = 1.0. 



that have a given momentum Pj as 



N,=L J2 IW(0)|A)P 



(11) 



In the LL model wc have Pj :— (27r/L)j. 

Solving the Bethe ansatz equations for a large number 
of eigcnstates we observe numerically that eigenvalues Nj 
arc given in decreasing order with respect to integer j: 
Nq > Ni > N2 > ■ ■ ■ ■ It thus follows that condensate 
fraction which corresponds to the largest eigenvalue of 
the one-particle reduced density matrix pi is indeed given 
by no = No/N, where A^o has been defined by sum (fTT|) 
over all eigcnstates with zero momentum. 

The estimates of condensate fraction no are plotted 
against coupling constant c in the upper panel of Fig. [1] 
over a wide range of c such as from c = 10~^ to c = 10'^ 
for different values of particle number N such as iV ~ 4, 
10, . . . , 400. For each N, condensate fraction uq becomes 
1.0 for small c such as c < 0.01, while it decreases with 
respect to c and approaches an asymptotic value in the 
large c region such as c > 100 or 1000. The asymptotic 
values depend on particle number N for N = 4, 10, . . . , 
400, and they are consistent with the numerical estimates 
of occupation numbers for the impenetrable ID Bose gas 
(see eq. (56) of Ref. 0). In the lower panel of Fig. [l] 
we plot fractions Uj for j = 0,1 and 2 against coupling 
constant c from c = 10"^ to c = 10^ with N = 20. The 
asymptotic values of rij for large c (i.e. c = 1000) are 
consistent with the numerical estimates for the impene- 
trable ID Bose gas (for ni and 712, see eqs. (57) and (58) 
of Ref. 0, respectively). 

We observe that condensate fraction no decreases as 
particle number TV increases when density n = N/L is 
fixed. It is the case for c < 0.1 in the upper panel of 
Fig. 1. Condensate fraction tiq decreases as TV increases 
even for small c such as c = 0.01, as shown in Fig. 2. 
Thus, it is necessary for coupling constant c to decrease 
with respect to A'^ so that condensate fraction no remains 
constant as TV increases with fixed density n. 

We now show the finite-size scaling of condensate frac- 
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FIG. 3: (Color online) Contours of condensate fraction no are 
plotted for various values of no in the 7 versus 1/N plane. 
Each contour is approximated by (|12|) : 7 as a function of 
1/N. 

tion riQ. In Fig. [3] each contour line gives the graph 
of interaction parameter 7 as a function of the inverse 
of particle number N for a fixed value of condensate 
fraction ng. They are plotted for various values of tlq 
from no = 0.6 to 0.99, and are obtained by solving the 
Bethe-ansatz equations numerically. For different values 
of density such as n = 1, 2 and 5, we plot countour lines 
with fixed values of condensate fraction no in the plane 
of interaction parameter 7 versus inverse particle num- 
ber We observe that the contours with the same 
condensate fraction no for the different densities coincide 
in the 7 versus 1 /N plane and are well approximated by 

7 = A/N"^. (12) 

Thus, condensate fraction no is constant as particle num- 
ber N becomes very large if interaction parameter 7 is 
given by the power of particle number N as in eq. (|12p . 

Applying the finite-size scaling arguments, we suggest 
from eq. (|12p that condensation fraction no is given by 
a scaling fuction </>(•) of a single variable ^N^: uq = 
(p^jN^). Here we recall the coincidence of contours for 
the different values of density n in Fig. [31 We thus 
observe that exponent rj and amplitude A of eq. 
are determimed only by condensate fraction no and are 
independent of density n. 

Let us consider amplitude ^ as a function of no. Wc 
denote it by ^ = /(no). Then, the scaling function (/){■) 
is given by the inverse function: no = f~^{A). In Fig. 21 



FIG. 4: (Color online) Exponent rj and amplitude A as func- 
tions of condensate fraction ?io ■ 

exponent rj increases with respect to no, and amplitude 
A decreases monotonically with respect to no. 

It follows from (fT2l) that BEC docs not occur in the 
ID Bose gas if wc fix parameter 7 and density n as sys- 
tem size L goes to infinity. However, if 7 is small enough 
so that it satisfies eq. (|12p for a given value of conden- 
sate fraction no, the ID Bose gas shows BEC from the 
viewpoint of the Penrose and Onsager criterion. Wc sug- 
gest that if condensate fraction no of a quantum state is 
nonzero and finite for large N, the mean-field approxima- 
tion is valid for the quantum state. For instance, there 
exist such quantum states that correspond to classical 
dark solitons of the GP equation Q, if parameter 7 is 
small enough so that it satisfies ([T2|) . 

With the scaling behavior wc derive various ways 
of the thermodynamic limit such that condensate fraction 
no is constant. For instance, we consider the case of a 
finite particle number, N = iVf. Choosing a value of no, 
we determine 7 by eq. (|12p as 7 = A{no)/Nf^"°\ Then, 
the ID Bose gas with N = Nf has the same condensate 
fraction no for any large value of L if coupling constant 
c is given by c = A{no)Nf~'^ / L. Let us set 77 = 1 and 
A'^f = 10, for simplicity. We have no — 0.97 in Fig. SI 
and 7 = 0.3 at 1/iV = 0.1 in the contour of no = 0.97 in 
Fig. m By assuming n = 1, it corresponds to the case 
of L = 10 and c = 0.3, and we have A = cL = 3, which 
is consistent with Fig. |4l Therefore, the ID Bose gas 
with N{ = 10 has no = 0.97 for any large L if c is given 
by c = 0.3/ L. Moreover, we may consider other types of 
thermodnamic limits. When density n is proportional to 
a power of L as L", condensate fraction no is constant 
as L goes to infinity if we set c oc 
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The scaling law (|12p and the estimates of condensate 
fraction in the present Letter should be useful for esti- 
mating conditions in experiments of trapped cold atomic 
gases in one dimension For instance, we suggest 

from Fig. 1 that BEC may appear in ID systems with a 
small number of bosons such as = 20 or 40 for c = 1 
or 10. 

In conclusion, we exactly calculated the condensate 
fraction of the ID Bose gas with repulsive interaction 
by the form factor expansion. We have shown the finite- 
size scaling behavior such that condensate fraction no is 
given by a scaling function of interaction parameter 7 
times some power of particle number N: uq = (pipfN^). 
Consequently, if parameter 7 decrease as 7 = A/N'^, con- 
densate fraction no remains nonzero and constant as par- 
ticle number N becomes very large. By modifying the 
thermodynamic limit, the ID Bose gas shows BEC from 
the viewpoint of the Penrose- Onsager criterion. 
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